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ABSTRACT 

We present the first measurements of the abundances of a-elements (Mg, Si, and S) extending 
out to beyond the virial radius of a cluster of galaxies. Our results, based on Suzaku Key Project 
observations of the Virgo Cluster, show that the chemical composition of the intra-cluster medium 
is consistent with being constant on large scales, with a flat distribution of the Si/Fe, S/Fe, and 
Mg/Fe ratios as a function of radius and azimuth out to 1.4 Mpc (1.3 r 2 oo)- Chemical enrichment 
of the intergalactic medium due solely to core collapse supernovae (SNcc) is excluded with very high 
significance; instead, the measured metal abundance ratios are generally consistent with the Solar 
value. The uniform metal abundance ratios observed today are likely the result of an early phase of 
enrichment and mixing, with both SNcc and type la supernovae (SN la) contributing to the metal 
budget during the period of peak star formation activity at redshifts of 2-3. We estimate the ratio 
between the number of SN la and the total number of supernovae enriching the intergalactic medium 
to be between 12-37%, broadly consistent with the metal abundance patterns in our own Galaxy or 
with the SN la contribution estimated for the cluster cores. 

Subject headings: ISM: abundances — X-rays: galaxies: clusters — galaxies: clusters: intracluster 
medium — galaxies: clusters: individual (Virgo Cluster) 


1. INTRODUCTION 

The chemical composition of the diffuse hot gas in 
galaxies and clusters of galaxies contains valuable infor¬ 
mation about the origin of chemical elements and their 
distribution during the evolution of the Universe. Ele¬ 
ments heavier than and including oxygen are produced 
by supernova explosions, with core-collapse supernovae 
(SNcc) producing predominantly lighter metals, from O 
to Si and S, while type la supernovae (SNIa) are r espon¬ 
sible for heavier e lements from Si to Fe and Ni (e.g. Tsuji-, 
moto et al.|1995 1. These elements are then expelled from 
their host galaxies by galactic winds and ram-pressure 
stripping. The metal abundances in the intergalactic 
space are therefore very sensitive to the numbers of SN la 
and SNcc that contributed to the enrichment, as well as 
to the time-scale over which the supernova products were 
injected into the intra-cluster medium (ICM). The initial 
metallicity of the progenitors, the initial mass function 
(IMF) of the stars that explode as SNcc, and details 
about the exact SNIa explosion mechani sm, also influ¬ 
ence the observed abundance ratios (see Werner et al. 


SN la products at the cluster centers (e.g. De Grand! fc 


2008 for a review). 

Initial results from the ASGA satellite suggested a radi¬ 
ally increasing Si/Fe ratio in the ICM, with the X-ray gas 
at radii beyond about 1 Mpc having a chemical composi¬ 
tion c onsistent with pure SNcc yields (Finoguenov et al. 
20001. The paradigm that emerged from these measure- 


ments was that of an early enrichment driven predom¬ 
inantly by SNcc in the proto-cluster phase, which led 
to the SNcc products being well mixed throughout the 
ICM, and a later continuous enrichment by SN la with 
longer delay times in the cD galaxy, creating a peak in 
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these measurements are concentrated within the central, 
brighter regions of any given cluster, therefore probing 
only a fraction of the total volume of the ICM. 

More recently, using a detailed s tudy of the outskirts 
of the Perseus Cluster with Suzaku, Werner et al. (20131 
presented measurements of the Fe abundance in the ICM 
in 78 independent regions beyond the cluster core and 
out to the virial radius (defined here as r 2 oo, the ra¬ 
dius within which the mean enclosed mass density of the 
cluster is 200 times the critical density of the universe 
at the cluster redshift). However, the equivalent widths 
of spectral lines from other elements are much weaker 
than Fe and therefore alpha-elemen ts have been s ignifi- 


cantly more challenging to detect. Sasaki et al. (20141 


reported nearly constant alpha-element to h'e ratios con¬ 
sistent with the Solar value out to 0.5 rigo in four nearby 
galaxy groups, but the metal abundance ratios beyond 
about a half of the virial radius, as well as the chemi¬ 
cal composition in the outskirts of more massive systems 
such as galaxy clusters, have not yet been measured. 

Here, we present results from Suzaku Key Project ob¬ 
servations of the Virgo Cl uster, the neyest g alaxy cluster 
at a distance of 16.1 Mpc (Tonry et al.|2001) and the sec- 
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ond brightest extended extragalactic soft X-ray source. 
Because line emission from Si and S is maximised around 
the mean temperature of the Virgo ICM 2 keV), this 
system presents an ideal target for determining the abun¬ 
dances of these metals. 

2. OBSERVATIONS AND DATA REDUCTION 

A large mosaic of Suzaku observations of the Virgo 
Cluster was obtained as a Suzaku Key Project during 
AO-7 and AO-8, totalling 60 pointings and 1 Ms of net 
exposure time. The observations cover the cluster o ut to 
beyond its estim ated virial radius (r 2 oo 1-1 Mpc, lUr- 
ban et al. 2011) along four arms extending towards the 
north, west, south, and east from the brightest cluster 
galaxy, M87. In addition, we included seven fields cover¬ 
ing the outer parts of the cluster along the northern arm, 
between 1.0 and 1.5 Mpc (AO-6). 

The data from the X-ray Imaging Spectrometers (XIS) 
0, 1 and 3 were analyzed fo ll owing the procedure de- 


scribe d in Simionescu et al. (2013) and |Urban et al. 
(2014). In short, in addition to applying the screening 
criteria recommended by the instrument tearrQ we fil¬ 
tered the data to include observation periods with the 
geomagnetic cut-off rigidity COR > 6 GV, and excluded 
two columns on either side of the charge-injected columns 
in XISl to avoid charge leak known to affect these areas 
of the detector. Furthermore, we used a mask file pro¬ 
vided by the XIS team to manually filter out pixels that 
are known to be hot or flickering We applied the latest 
contamination layer calibration from 2014 August 25. 

An image of the complete mosaic in the 0.7-2.0 keV 
energy band, corrected for instrumental background and 
vignetting, is shown in Figure[2 Point sources were iden¬ 
tified and excluded from further analysis using the CIAO 
tool wavdetect employing a single wavelet radius of 1 ar- 
cmin, which is matched to the Suzaku half-power radius. 
In addition, we removed a 13' radius circle around the 
background cluster Abell 1553. 

3. SPECTRAL MODELING 

For each arm, we extracted spectra from partial annuli 
centered on the core of M87. We imposed a minimum 
of 1200 source counts per spectral extraction region, ac¬ 
counting for the i nstr umental and sky backgrounds esti¬ 
mated in Section |3.1[ Due to the point-spread function 
(PSF) of Suzaku, the minimum width of each annulus 
was set to 3 arcmin. Along the southern arm, we only 
use the data out to a radius of 957 kpc, corresponding 
to the minimum surface brightness between M87 and the 
infalling halo of M49. 

The spectral modelling was performed with the 


XSPE C spectral fitting package (version 12.8.2, Arnaud 
1996), employing the modified C-statistic estimator. The 


spectra were binned to a minimum of one count per chan¬ 
nel and fitted in the 0.7-7.0 keV band. Unless noted 
otherwise, the ICM was modelled as a thermal plasma 
in collisional ioniza tion equilibrium using the vapec code 
(Smith et al.||200lj. The metal abundances are reported 


^ Arida, M., XIS Data Analysis, 
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Fig. 1. — Image of the complete Suzaku mosaic of the Virgo 
Cluster in the 0.7-2.0 keV energy band, corrected for instrumental 
background and vignetting. The white dashed circle marks the 
approximate location of r2oo? while the red circles mark the RASS 
regions used to estimate the Galactic foreground flux. The region 
excluded around the background cluster Abell 1553 is shown in 
yellow. Colorbar units are counts/s/arcmin^. 


with respect to the proto-solar units of Fodders (2003). 
With res pect to more recent S olar abundance measure- 
ments by|Asplund et al.| (|2005|), the Si/Fe would be 14% 
higher, S/h'e higher, and Mg/Fe 5% lower than the 
values reported here. The Galactic absorption column 
density, uh, was fixed to the average value at the loca¬ 
tion of each respective Suzaku fie ld based on the radio 
HI survey of Kalberla et al. (2005). 


3.1. Background subtraction 

To account for the instrumental and cosmic back¬ 
ground we f ollow the methods laid out in e.g. |Urban| 
(2014). In short, the instrumental background 

way, using the task 


et al. 


was 


In short, 
subtracted in the standard 


xisnxbgen (Tawa et al. 2008). The cosmic X-ray back¬ 
ground model consisted of a power-law emission compo¬ 
nent that accounts for the unresolved population of point 
sources (CXB), one absorbed thermal plasma model for 
the Galactic halo (GH) emission, and an unabsorbed 
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TABLE 1 

Background model parameters. The first set of errors 

GIVES statistical UNCERTAINTIES AT 68% CONFIDENCE; THE 
SECOND SET REPRESENTS SYSTEMATIC UNCERTAINTIES DUE TO 
SPATIAL VARIATIONS. 



kT (keV) 
or r 

arm 

XSPEC norm 
per arcmin^ (xl0“^) 

CXB 

1.53 ± 0.04 

all 

0.83 ±0.05 ±0.06 



N 

1.31 ±0.10 ±0.22 

GH 

0.20 ± 0.01 

W 

1.14±0.10±0.12 



s 

2.20 ±0.15 ±0.66 



E 

1.73 ±0.12 ±0.11 



N 

1.08 ±0.04 ±0.08 

LHB 

0.104 ± 0.002 

W 

1.35 ±0.05 ±0.05 



s 

1.46 ±0.06 ±0.27 



E 

1.41 ±0.05 ±0.11 


thermal plasma model for the Lo cal Hot Bubble (LHB) 
(see e.g. Kuntz fc Snowden||2000 ). 

The average hux and spatial variation of the GH and 
LHB emission was determined using a set of 12 ROSAT 
All-Sky Survey (RASS) spectra obtained from circular 
regions with a radius of 1 degree (three regions beyond 
each Suzaku arm; see Figure H. The metallicity of both 
components was assumed to be Solar. Six Suzaku point¬ 
ings (the outermost three along the N and W directions, 
respectively) do not contain any significant ICM emis¬ 
sion. We thus determined the parameters of the CXB 
model by fitting the spectra obtained from the full field 
of view of these observations in the 2-7 keV band. The 
best-fit model parameters and their estimated statistical 
and systematic uncertainties are summarised in Tabled 

We checked for solar wind charge-exchange (S WCX) 
emission contamin ation following the analysis of |Fuji-| 


moto et al. (2007). Based on solar wind proton flux 
measurements from the WIND SWF (Solar Wind Ex- 


perimenl|^, in combination with estimates of the radius 
of the Earth’s magneto-sheath along the Suzaku line of 
sight during our observations, we find that the level of 
SWCX contamination is expected to be negligible for all 
the pointings considered in the present analysis. 


4. RESULTS 

In this Letter, we focus on measurements of the abun¬ 
dance of alpha-elements in the outskirts of the Virgo 
Cluster. A more detailed discussion regarding the ther¬ 
modynamical properties of the ICM, and the total gas 
and metal masses compared to the stellar component, is 
the subject of a separate publication (Simionescu et ah, 
in prep). 

We divided the data into four radial bins, spanning 
150-200 kpc, 200-300 kpc, 300-500 kpc, and from 500 
kpc out to the edge of the detection limit. For each radial 
bin, we fit in parallel the spectra from all partial annuli at 
any azimuth that fall within the respective radial range. 
The temperature and spectrum normalization were left 
as free parameters for each partial annulus, while the 
Fe, Si, S, and Mg abundances were coupled to the same 
value. The abundances of all other elements between O 
and Ni were set to 0.2 Solar. The results are summarised 
in Tabled 


^ http://web.mit.edu/afs/athena/org/s/space/www/wind.html 


TABLE 2 

Metal abundances as a function of radius for various 

SPECTRAL MODELS. ERRORS ARE QUOTED AT AC = 1. 



150-200 kpc 

200-300 kpc 

300-500 kpc 

500-1400 kpc 



vapec 

IT 


Mg 

0.16 ±0.06 

0.33 ±0.09 

0.15 ± 0.09 

0.14 ± 0.09 

Si 

0.14 ± 0.03 

0.17 ± 0.05 

0.15 ± 0.05 

0.10 ±0.05 

S 

0.21 ± 0.05 

0.22 ±0.07 

0.28 ±0.08 

0.21 ± 0.08 

Fe 

0.22 ±0.01 

0.28 ±0.02 

0.15 ± 0.02 

0.17 ± 0.02 



vapec 

2T 


Mg 

0.27 ± 0.07 

0.50 ±0.11 

0.39 ±0.12 

0.32 ±0.11 

Si 

0.18 ± 0.04 

0.24 ± 0.05 

0.23 ±0.07 

0.15 ±0.06 

S 

0.24 ± 0.05 

0.28 ±0.08 

0.33 ±0.10 

0.25 ±0.10 

Fe 

0.24 ± 0.01 

0.32 ±0.02 

0.28 ±0.02 

0.23 ±0.03 



spex 

2T 


Mg 

0.21 ± 0.07 

0.48 ±0.10 

0.32 ± 0.12 

0.45 ±0.12 

Si 

0.20 ± 0.04 

0.29 ±0.06 

0.25 ± 0.05 

0.21 ± 0.07 

S 

0.24 ± 0.05 

0.29 ±0.08 

0.33 ±0.11 

0.21 ± 0.12 

Fe 

0.29 ±0.02 

0.39 ±0.03 

0.32 ± 0.04 

0.32 ±0.03 


Biases in measuring the metallicities in the ICM, par¬ 
ticularly that of Fe, can arise when X-ray emitting gas 
at several different temperatures contributes to the in¬ 
tegrated spectrum, especia.lly if the average temperature 
of the ICM is low (Buote 2000). We have thus refitted 
the data with a two-temperature model (hereafter 2T 
model or 2T fit), with the abundances of the two gas 
phases assumed to be the same, the normalisation of the 
cool component left free in the fit for each spectral ex¬ 
traction region separately, and the temperature of the 
cooler gas, kTcooi, successively fixed at 0.6, 0.8, 0.9, 1.0, 
1.1, and 1.2 keV. The largest differences compared to 
the IT fit are obtained for kTcooi=^ keV, and we show 
the results for this model in Table [21 All metal abun¬ 
dances are systematically higher than those derived from 
a single-temperature ht; however, the metal abundance 
ratios are more robust, with Si/Fe and S/Fe largely con¬ 
sistent between the IT and 2T models, while Mg/Fe is 
systematically higher by approximately Icr for the two 
temperature fit in each radial bin compared to the single 
phase approximation. 

Given that it is likely that multi-temperature struc¬ 
ture is present in the ICM, and that more complex multi¬ 
phase models can not be constrained (and are unlikely to 
have a bigger impact than the difference between single¬ 
phase and 2T models), we have chosen to use the re¬ 
sults from the 2T fit as reference values for later discus¬ 
sion. For this spectral model, in the radial bin beyond 
500 kpc, we detect Mg with a statistical significance of 
2.9(7, and Si and S each with a significance of 2.5cr; line 
emission from alpha-elements in the cluster outskirts is 
thus detected with a cumulative significance of > 4.5 (t. 
We have verified that the same level of significance is 
obtained when the statistical uncertainties are derived 
using a Markov chain Monte Carlo approach instead of 
the AC = 1 method adopted throughout the rest of this 
manuscript. 

In the left panel of Fig. we show the metal abun¬ 
dance ratio profiles obtain^ from the 2T fit. All the 
alpha-element to Fe ratios are consistent with being con¬ 
stant as a function of radius all the way to the edge of the 
Virgo Cluster. The best-fit constants in the radial range 
between 150 and 1400 kpc are 0.75±0.10, 0.99±0.14, and 
1.31 ±0.19 for the Si/Fe, S/Fe and Mg/Fe ratios, respec- 
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Fig. 2.— Left: Abundance ratio pr ofiles from the 2T vapec fit with kTcooi = 1-0 keV are shown in blue. Red data points show published 
results from|Simionescu et al.|(|2010|l for the N arm. Systematic uncertainties due to background fluctuations are shown with dotted lines 
for the GH and dashed lines for the CXB power-law. Right: Abundance ratios as a function of azimuth (measured counterclockwise of 
west). Dotted horizontal black lines show the Icr confidence interval for the best-fit constant to the corresponding radial profile. Magenta 
dashed lines show the expected values for enrichment by pure SNcc yields in both plots. 


tively, therefore generally in good agreement with the 
Solar abundance pattern, with slightly sub-solar Si/Fe 
and super-solar Mg/Fe. Linear hts to the radial profiles 
of the abundance ratios yield a slope consistent with zero 
in all cases. Our results are in good agreement with the 
metallicity ratios measure d within the centra l O.drigo for 
four other galaxy groups ( Sasaki et al.||2014 ). 

We have also tested for the azimuthal dependence of 
the abundance ratios by fitting in parallel the spectra 
from all radii along a given Suzaku arm. The results 
are shown in the right panel of Figure We find that 
the measured abundance ratios are consistent with a flat 
distribution as a function of azimuth; along each arm, the 
measurements agree with the best-ht constant obtained 
from the corresponding radial prohles. 

4.1. Systematic uncertainties 

To investigate the systematic uncertainties associated 
with the assumed plasma emission code, we have also re¬ 
fitted the data using the CIE p lasma model implemented 
in spex (Kaastra et al. 1996b version 2.00.11. The re¬ 


sults for a two-temperature fit with kTcooi fixed at 1.0 
keV are shown in Table |2] The abundances of all the 
alpha-elements agree with those derived from the vapec 
2T fit within the Icr statistical error bars; however, due 
to differences in the Fe L-shell transitions between the 
two plasma codes, the Fe abundance is systematically 
higher when using spex. This only increases the signif¬ 
icance with which a contribution from SN la products 
would be required to explain the enrichment of the clus¬ 
ter outskirts (Section]^. 

We estimated the systematic uncertainties due to the 
background model by refitting the spectra with the GH 
and CXB fluxes scaled according to the expected fluc¬ 
tuations described in Table accounting for the I/VH 


dependence of the cosmic variance of unresolved point 
sources on the size of the extraction region H. The re¬ 
sults are shown in Figure!^ The systematic uncertainties 
are in all cases smaller than the statistical error bars of 
the measurements. The contribution of the LHB and the 
associated systematic uncertainties are negligible in the 
energy band considered here (above 0.7 keV). The effects 
of a 5% uncertainty on the instrumental background are 
comparable to those due to CXB fluctuations. Spatial 
variations of the Galactic column density do not influ¬ 
ence the results significantly, such that consistent results 
are obtained assuming a uniform uh across the entire 
Suzaku mosaic. 

Note that, although the conclusions based on the Si 
abundance measurements provide a high statistical con¬ 
fidence, they may carry additional systematic uncertain¬ 
ties compared to S and Mg due to the complexity in cal¬ 
ibrating the Si instrumental line as well as the effective 
area around the Si edge. 

5. DISCUSSION AND CONCLUSIONS 

We present the first measurements of the abundances 
of a-elements extending out to beyond the virial radius 
of a galaxy cluster, specifically the nearby Virgo Gluster. 
Our main findings are summarised below. 


5.1. Pure SNcc enrichment ruled out in the cluster 
outskirts 


Assuming the SNcc yields of Nomoto et al.|(2006) for 
Solar metallicity of the progenitor and averaged over a 
Salpeter IMF, the expected Si/Fe, S/Fe and Mg/Fe ratios 
from enrichment of the ICM due solely to core-collapse 
supernova products are 2.31, 1.92, and 2.93 Solar, re¬ 
spectively (shown as dashed magenta lines in Fig[^ these 
values increase if the metallicity of the progenitor is sub- 
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solar). Therefore, the best-fit constant values describing 
our measured alpha-element to Fe abundance ratio pro¬ 
files rule out pure SNcc enrichment beyond the cluster 
core with very high significance (15.6cr, 6 . 6 cr, and 8.5a 
for Si/Fe of 0.75 ± 0.10, S/Fe of 0.99 ± 0.14, and Mg/Fe 
of 1.31 ±0.19, respectively). Using only the data beyond 
half of the cluster’s virial radius, pure SNcc enrichment 
is excluded at the 5.8a level based on the best-fit Si/Fe 
ratio, and l.Scr and 2.7cr for the S/Fe and Mg/Fe mea¬ 
surements, respectively. 

5.2. Uniform chemical composition of the Universe over 
a wide range of scales 

Rather than following the chemical enrichment pat¬ 
terns expected from pure SNcc enrichment, as inferred 
from earlier work, the measured metal abundance ratios 
are largely consistent with the Solar value over the en¬ 
tire cluster volume probed by the present measurements 
(both as a function of radius and azimuth, from just out- 
side the influence o f M87 out to 1 . 3 r 2 oo)- In addition, 
Million et al. (201l|) report detailed metal abundance ra- 
tios in the M87 halo using deep Chandra data, and find 
a Mg/Fe ratio also consistent with the Solar value from 
a radius of 7-40 kpc, while the Si/Fe and S/Fe ratios 
decrease from 1.2 Solar at a radius of 10 kpc to 1.0 Solar 
at 40 kpc. Generally, therefore, given the wide range of 
spatial scales probed, from the chemical content of stars 
like our Sun to the enrichment pattern on megaparsec 
scales, the chemical composition of the Universe appears 
to be remarkably uniform. 

This supports the conclusions based on previous 
Suzaku Key Project observations of the Perseus Cluster, 
which revealed a uniform Fe abundance distribution on 
large scales, indicating that most of the metal enrichment 
of the intergalactic medium occurred before the entropy 
in the ICM became stratified, probably during the pe¬ 
riod of maximal star formation and blac k hole activity 
at redshifts z = 2 —3 ( Werner et al.|20T3 ). Such an early 
enrichment and mixing is likely responsible for producing 
the uniform metal abundance ratios observed today. 

The average Fe abundance in the outer parts of the 
Virgo Cluster is 0.23 ±0.03 Solar, marginally lower com- 
pare d to 0.29 ± 0.01 in the Perseus Cluster (converted 
from Werner et al.|2013 to the Solar units adopted here). 
This ditterence may be due to uncertainties in the Fe- 
L line emissivities, since the spex plasma emission code 
suggests a higher Fe abundance that is in agreement with 
the Perseus Cluster measurement. 


5.3. The relative contribution of SN la to the early 
enrichment of the intergalactic medium 


In order to provide a homogeneous chemical composi¬ 
tion of the Universe, the rapid production of a-elements 
by SNcc associated with the peak star formation activity 
at redshifts of 2-3 must have been complemented by a 
sufficient number of SN la exploding during or shortly 
after this epoch. This is consistent with results from 
type la supernova delay-time distributions, which imply 
that a large fraction of such supernovae explode less than 
^ 0.5 Gyr after the formation of the progenitor binary 
system (Mannucci et al.||2006|. 

The ratio between the number of SN la and the total 
number of supernovae contributing to the enrichment of 
the ICM, R%ia = Nia/{Nia ± Ncc): can be estimated 
by mixing model yields of SN la and SNcc products in 
various percentages until the observed metal abundance 
ratios are reproduced. Using the same SNcc yields men¬ 
tioned above, an d assuming the SN la y ields based on the 
WDD2 model of |Iwamoto et al.| ( |1999[ ), we obtain R%ia 
values of 37tg^%, 181 / 5 %, ^-’^cl 12 + 2 ^, corresponding to 
the measured Si/Fe, S/Fe, and Mg/Fe, respectively. Note 
however the additional possible systematic uncer taint ies 
associated with the Si/Fe determination (Section |4.1[ ). 

Although additional uncertainties due to different su¬ 
pernova yield models and assumptions about the IMF 
and initial metallicity of the progenitor stars re main (see 
detai led discussion in e.g. |de Plaa et al.|[2007l |de Plaa| 
20131, the ratios presented here are generally con sistent 


with those estim ated for our own Galaxy (15%, Tsuji- 
moto et al.|1995[) or for the cluster cores (about 20-30% 


ai.||2007| ISimionescu et al.||2009H2010|). 

that the fractional contributions calc 


e.g. de Plaa et 


Finally, note tJiat tlie tractional contributions calcu¬ 
lated here only reflect the ratio of SN la to SNcc that 
enrich the ICM; if the process of injecting metals from the 
progenitors into the intergalactic medium has a different 
efficiency or timescale for different supernova types, the 
true R%ia may be different than reflected by the compo¬ 
sition of the ICM. 
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